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AN EXTENSION TO LOWER PRESSURES OF THE 


NAVTER-STOKES THEORY FOR THE VISCOUS FLOW OF GASES 
BETWEEN RELATIVELY ROTATING COAXIAL CYLINDERS 

By Willard E. Meador 
Langley Research Center 

SUMMARY 

An extension of the Grad 13-moment approximation is applied to the velocity- 
linearized Boltzmann kinetic equation for the problem of viscous gases flowing between 
two relatively rotating coaxial cylinders. Precise boundary conditions are employed in 
deriving a velocity distribution function which is exact through quadratic terms in the 
Knudsen number and which thus provides reliable generalizations of previous work to 
lower pressures; in particular, the Navier-Stokes expression for the torque experienced 
by the inner cylinder is shown to suffice for Knudsen numbers extending well into the 
transition flow regime if the separation between the cylinders is small compared with 
their radii. Low-pressure effects for larger ratios of these distances differ significantly, 
however, from the calculations of Lin and Street based on the Burnett equations. At 
more moderate pressures corresponding to the retention of terms through the square of 
the Knudsen number in the torque and through the first power of the Knudsen number in 
the velocity distribution function, the simple Navier-Stokes theory is always valid and 
consistent in the problem described. 

Results of the present research are compared also with those of Wang Chang and 
Uhlenbeck, who use the obsolete Maxwellian separation of reflected particles into specular 
and diffuse groups, and with previous formulations based on Chapman-Enskog expansions 
and the Grad 13-moment approximation. The more important conclusions are summa- 
rized as follows: (1) the velocity jumps at the cylinder surfaces generally follow 

Maxwell's original predictions rather than Wang Chang and Uhlenbeck's later modifica- 
tions, (2) the basic character of the velocity profile in the midst of the gas does not 
change as the surfaces are approached, and (3) 13 moments are insufficient to describe 
low-pressure systems. 


INTRODUCTION 

The problem of viscous gases flowing between two relatively rotating coaxial 
cylinders has received much attention both on its own merit and because such systems 


are especially suited for the measurement of velocity slip at solid surfaces. (See ref. 1, 
pp. 291-300, for the earliest work.) Despite the interest, however, several important 
questions remain unanswered concerning the precise kinetic behavior of a gas near a 
boundary and the range of validity of the Navier -Stokes expression (ref. 1) for the torque 
experienced by the stationary inner cylinder. In particular, the assumption that the 
velocity distribution function calculated in the midst of the gas extends uniformly to the 
cylinder surfaces needs to be investigated, as does the appearance of terms proportional 
to the square of the mean free path in macroscopic relations derived from the second 
Chapman-Enskog approximation (ref. 2). More recent research on similar subjects in 
the case of plane Couette flow is summarized in reference 3. 

A prime objective of the present research is to resolve these difficulties by exactly 
solving the corresponding velocity -linearized Boltzmann kinetic equation for Maxwellian 
particles through quadratic terms in the Knudsen number and imposing the conservation 
of particles, momentum, and energy at the boundaries. These conservation conditions 
replace the obsolete Maxwellian hypothesis (refs. 4 and 5) that separates the reflected 
particles into specular and diffuse groups and adds contributions to the velocity distribu- 
tion function which are nonanalytic in the mean free path (ref. 6). As explained in refer- 
ence 3, the Maxwellian hypothesis is nonphysical because it imposes the same accommo- 
dation or absorption coefficient for every macroscopic moment of the gas. The restric- 
tion to gases of Maxwellian particles (i.e., particles which attract or repel each other 
with forces varying as the inverse fifth power of their separation) is not expected to alter 
the basic conclusions. 


SYMBOLS 

~c particle velocity 

d intercylinder separation distance 

D drag force per unit area in limit of infinite cylinder radii 

f velocity distribution function 

f Maxwellian distribution function relative to local flow velocity 

f Q ^r^ Maxwellian distribution function relative to stationary inner cylinder 

f o ( r 2 ) Maxwellian distribution function relative to rotating outer cylinder 
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p 

~p 

. 0 . 

p 

Q 

r,e,z 


r l’ r 2 

,e 

t 

T 


u 


contributions to perturbation function 
arbitrary velocity function 

perturbation function relevant to gas at surface of inner cylinder 
perturbation function relevant to gas at surface of outer cylinder 
Boltzmann's constant 
mean free path 

torque per unit length of stationary inner cylinder 

particle mass 

particle number density 

scalar gas pressure 

pressure tensor 

traceless pressure tensor 

heat-flux tensor 

cylindrical coordinates (r measured from axis of coaxial cylinders); also 
signify vector and tensor components when used as subscripts 

radii of inner and outer cylinders, respectively 

unit vectors associated with cylindrical coordinate system 

time 


temperature 


particle velocity relative to frame of reference moving with local gas flow, 


nondimensionalized, 
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U unit tensor 

v local gas flow velocity 

Vj magnitude of gas flow velocity at surface of stationary inner cylinder; also 

called slip velocity or velocity jump 

Vg magnitude of gas flow velocity at surface of outer cylinder 

v,„ rotational speed of outer cylinder 

w 

v mean equlibrium particle speed 

X body force per unit mass 

accommodation coefficient (inner cylinder) for energy absorption 
y particle velocity relative to stationary inner cylinder, nondimensionalized, 



£ slip distance 

slip distances relative to inner and outer cylinders, respectively 

7 ] viscosity 

ju particle velocity relative to rotating outer cylinder, nondimensionalized, 



p mass density 

ct 1’°2 accommodation coefficients (inner and outer cylinders, respectively) for 

parallel momentum absorption 

accommodation coefficient (inner cylinder) for perpendicular momentum 
absorption 
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a'i accommodation coefficient (inner cylinder) for parallel heat-flux absorption 

(see eq. (34)) 

4> first-order perturbation function 

a) angular frequency of outer cylinder 

Special notation: 

collisional time derivative 



<G> 

< G | h l> 


velocity average of G referred to complete distribution function 

< G h> velocity averages of G referred to and f^rg^hg, 

respectively 


Vector and tensor quantities without arrows signify magnitudes or components. 

GENERAL CONSIDERATIONS 


The torque exerted on unit length of a stationary inner cylinder (radius r j, 
infinite length) by the shearing action of a viscous gas flowing between it and a rotating 
outer coaxial cylinder (radius infinite length) is given by 

L = - 2 "l 2 p r 9 ( r l) (1) 

where P ^ is the rd - component of the pressure tensor 

= nm ((c - v)(^ - v)) (2) 

Hence, the determination of P r ^ from kinetic theory is the essential task and comprises 
the bulk of the present research. 

Calculations of such quantities usually begin with the general time -independent 
kinetic equation (ref. 2) 


c • Vf + X 


_8f 

8c" 



( 3 ) 
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the mc-moment of which yields the macroscopic equation of motion 


V • + pU + nmvvj = 0 

_o 

in the absence of body forces. The traceless pressure tensor P is defined by 
o _ 

P = P - pU 


( 4 ) 


( 5 ) 


If no temperature gradients are applied between the cylinders, the scalar pres- 
sure p depends on the radial distance r only through effects induced by the friction- 
generated heat. Since such effects are proportional to the square of the Mach number, 

as are the diagonal elements of P in the midst of the gas, the direct expansion of equa- 
tion (4) gives 


o 

V • P = V • 



( 6 ) 


through first-order (i.e., linear) terms in the flow velocity. A reasonable tentative 
assumption (to be confirmed or rejected by the solution of the kinetic equation subject to 
boundary conditions) is that equation (6) applies also in the immediate proximity of the 
cylinder surfaces. 

Although equation (6) is an important differential relation, its utilization in the 
present problem is more indirect than direct because the integration constant P^ (n)> 
rather than the variation of P r g, with r, is the significant parameter in equation (1). 
What equation (6) does provide is a differential relation for the flow velocity v if P r g 
can be expressed in terms of v through the solution of the kinetic equation. The inte- 
gration constants in the solution of the differential equation for v must be related to 
boundary conditions at the surfaces and also to the relative rotation rate of the two 
cylinders; hence, the procedure for writing L in equation (1) as a function of measur- 
able quantities is outlined by this discussion. 

Unfortunately, however, the difficulties inherent in cylindrical geometries preclude 
exact closed-form first-order velocity distribution functions of the type found for simple 
Couette flow between parallel plates (ref. 3); as a result, additional precautions are nec- 
essary to ensure complete first-order descriptions of phenomena relevant to cylindrical 
problems. Other difficulties arise from the nonlinear character of the velocity profile, 
the spatial variations of which are a priori unknown and must be derived according to the 
procedure in the preceding paragraph. 
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KINETIC THEORY 


A first-order solution of equation (3) can be written in the convenient form 

f - f 0 a + 0) = f 0 [ + § p re u r"e + « ( r -“)] < 7 > 

where g is an unknown function, f Q is the Maxwellian distribution function 



and u is a dimensionless particle velocity defined by 



All thermodynamic quantities are constant through first order in the absence of applied 
gradients, and the flow velocity satisfies 

v = 0v(r) (10) 

in the present problem. 

Since 



the substitution into equation (3) of equations (7), (11), and (12) yields 
c • Vf » c • Vf Q + f c • V(p 
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with the aid of equation (6) and the well-known properties of collision integrals for 
Maxwellian particles (ref. 2). The viscosity r] correctly represents the detailed colli- 
sion dynamics. 

Equation (13) is called the velocity-linearized Boltzmann equation because quad- 
ratic and higher powers of v are neglected in its derivation. The complexities intro- 
duced by the cylindrical geometry are obvious: As r approaches infinity, corresponding 
to flat plates, the exact first-order solution is equation (7) with g set to zero. In the 
present problem, however, g is more complicated and must contain, at the very least, 
some seldom -used elements of the general heat-flux tensor 



More specifically, g can be expressed as 



in order to balance the ^u r ^ - iu^ju^ term in equation (13). 

The insertion of equation (15) into the collision integral of equation (13) gives the 
expression 




3P r 0i 

f2kT 

pr 

( m 


2 1„ 2\ l/2kT 

U r ' 3 U 0 h + 2 lT 


u • Vg 




+ 3Q 


rr0 




u 0 + 


i Kto) 



Accordingly, the P r0 term in <p is proportional to 


if dv _ v\ 
v \dr r ) 


and the Q rr g term is proportional to 

dv _ v\ 
rv \dr r) 


(16) 
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if the mean free path l and the mean equilibrium particle speed v are defined by 


l = 


2r7 

Pv 


( 17 ) 


and 



(18) 


This analysis permits a definite statement about the physical consequences of 
dropping the u • Vg contribution to equation (16) and the collisional time derivative 
of f Q gj to obtain 


/ dv 
\dr 


v \ 3P 

7j u r u 0- pr 


^(¥)% 2 - 
Pr 9 “r^Q rre (^) 1/2 (u r 2 -|V) 


u. 


(19) 


The expansion is twofold in the sense that terms proportional to 



or 


; 3 »» 

r i 2 ( r 2 - r i) ' 


are neglected, as are higher powers of l and the Mach number v,.,/v of the outer 

w/ 

cylinder at r = rg. (Note that this definition of the Mach number differs by a constant 
factor from the one used in aerodynamics.) From this point forward the inner cylinder 
is considered to be stationary. 

Equations (7), (15), and (17) to (19) combine to give the following first-order results, 
which are complete through terms containing the square of the mean free path: 


and 



0 = 


2P 


re 


IT 




( 20 ) 


( 21 ) 


( 22 ) 
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These expressions, together with the differential relation 

d^v | 1 dv v _ Q 
<^ r 2 r dr r 2 


(23) 


obtained from equations (6) and (20), complete the kinetic-theory description of the gas 
between two relatively rotating coaxial cylinders. They extend the Navier-Stokes repre- 
sentation to larger Knudsen numbers through the appearance of terms in the velocity 
distribution function; in addition, they are derived from a kinetic equation incorporating 
the divergence of the traceless pressure tensor on the left-hand side, which is not normal 
procedure in the second Chapman-Enskog approximation (ref. 2). 

In line with the discussion following equation (6), the assumption is now made that 
equations (20) to (23) are valid in all regions occupied by the gas, including the close 
proximity to the cylinder surfaces. The only physical criterion for this assumption is 
whether the boundary conditions can be satisfied: If so, the description is accurate 
through first order in the flow velocity and through squares of the Knudsen number; if 
not, equation (22) is inadequate and higher moments must be considered. 


BOUNDARY CONDITIONS 


A convenient formulation of equations (7) and (22) for the application of boundary 
conditions at the surface /r = r^J of the stationary inner cylinder is obtained from the 
following manipulation of equation (8): 


° y2ffkT / 


m_ c 2 


2kT 


2vCg + 


m 


\3/2 


27rkTj 


e-r' 


l + 2[JL\ 1/2 vy, 


2kT J 


' f o( r l) f 1 


4v 

+ 7 — 

1 / 2 - 9 
7T V 


(24) 


where 



(25) 


Accordingly, the first-order velocity distribution function becomes 


f = f o( r l)( 1 + h l) 


(26) 
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with 


h l = 


4v 


2P 


n^ 2 v 


y 0 + 


rd 


y + 
' r 


_i* (V 2 -ly 2 

l/2 r V r 3 /0 


y 0 


(27) 


The introduction of h 1 permits the conservation laws for particles, the r-, 9-, 
and z-components of momentum, and energy to be expressed at r = r ^ by the relations 


°= <V r | 


(28) 


»N 1 <’' r 2 h> vo -<’'r 2 K> rr<0 -<V r 2 K) yr: 

<'i<V'e| h i> <0 = <V' 9 | h i> 

r r 

<’l<’V'z| h l> <0 ‘ (V'zh) 

r r 

and 

“l( y 2 y r| h l) <0 = ( y 2 y r| h l) 
y r 


(29) 


(30) 


(31) 


(32) 


respectively, where 0, Up and are accommodation coefficients and 

the "bra-ket" notation refers to velocity averages symbolized by the definition 


( G ( r l’ y )| h l) 4K(ri) h i G (^ ( 33 ) 

Each of equations (28) to (32) is interpreted in the same manner, the left-hand 
sides representing the average amounts of the corresponding macroscopic properties 
absorbed by unit surface area in 1 second and the right-hand sides representing the net 
fluxes toward the surface of the same properties. As explained in reference 3, these 
conditions replace the obsolete Maxwellian hypothesis that the reflected particles can be 
separated into diffuse and specular groups relating to a single accommodation coefficient 
for all moments. The coefficients are generally different for different moments. 
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Also mentioned in reference 3 is the fact that additional steady-state conservation 
laws are superfluous in the sense that they do not constrain the distribution function 
beyond what is implicit in equations (28) to (32) . For example, the result of imposing 
the he at -flux condition 


"i( r 2 r r r e| l 'i) j , <0 = <’'V«| h i) 


(34) 


on equation (27) is merely to write in terms of o^, which does not affect the solution 
of the kinetic equation in any way because equations (28) to (32) completely define (on the 
average) the gas-surface interactions. Yet a restriction similar to equation (34), but 
with cr'j incorrectly equated to according to Maxwell's hypothesis, was imposed by 
Wang Chang and Uhlenbeck (ref. 6) and led to complicated nonanalytic (in the mean free 
path) contributions to both the distribution function and the calculated macroscopic prop- 
erties. Such contributions are not implied by the more general and more correct bound- 
ary conditions employed in the present research. 

Of the five constraints listed in equations (28) to (32), only equation (30) is satisfied 
nontrivially by the perturbation function in equation (27). This condition provides a rela- 
tion between the flow velocity v^ at r = r ^ and the accommodation coefficient 
according to the following development: 


(?r y 8 


lh \ = _A 

V 1 

frPre^l) 

(i - iL\ 

1 V y r <0 " 

V 

4p 

l ”4 


<V'e| h i) 



(35) 


(36) 


and 


V 1 _ 7r ^l P r0 1 

( r i), 

( 1+ \ 

v 4p l \ 

V ^l r l/ 


where satisfies 

, ( 2 - ff i) z 

n ol 


(37) 


(38) 


Note that to the slip velocity v^ in equation (37) is added a new contribution (the 
second term) which originates from the heat-flux tensor element Q rr g (instead of the 
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vector heat flux which vanishes according to the u^u-moment of eq. (22)) and thus is not 
predicted by the Grad 13-moment velocity distribution function. Because of the explicit 
cancellation of the slip distance ^ and except for the implicit dependence of P^r-^ 
on the same parameter, this term enhances the slip velocity independently of the value 
of the accommodation coefficient cr^; in particular, the slip velocity is nonvanishing even 
when the average momentum component parallel to the cylinder surface is exactly 
reversed = 2) by the reflection. Although the precise physical explanation of this 
phenomenon is not clear, it can be called a geometric effect with some justification 
because of the r^ in the denominator. As r^ approaches infinity, corresponding to a 
transition from cylindrical to rectilinear gas flow, the effect becomes negligibly small. 

The only previous derivations of the geometric effect, which was not singled out or 
discussed as such, were those of Schamberg (ref. 7) and Lin and Street (ref. 8) on the 
basis of the third Chapman-Enskog approximation corresponding to the Burnett equations. 
Both references give slip velocities which differ from equation (37) by the substitution of 
the factor 10/3 for the factor 4 in the second term within the parentheses. Of more 
importance, however, is the fact that references 7 and 8 predict the form of the first- 
order expression for the velocity slip to be independent of the interparticle interaction 
potential; hence, except for the numerical value of the viscosity coefficient, which is 
somewhat irrelevant to the present research, strong support is provided for the assertion 
that the formal results of this paper are not restricted by the limitation of the detailed 
development to Maxwellian molecules. 

A similar study relating to the surface ^r = rgj of the outer cylinder begins with 
the first-order distribution function given by the expressions 


f = f. 


»( r 2)( 


1 + h 2) 


f o( r 2 ) n ( 27rk T) 


3/2 2 

e-ti 


and 


‘( v - *w) + 

IT 1/2 V 8 


2P 


rQ 


+ 


4 1 


-l/2-l'~ r 


a 2 _ A„ 2 

Mv 3 M 0 




where v w is the rotational speed of the surface and jn satisfies 



(39) 

(40) 


( 41 ) 


(42) 
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and continues with the application of the conservation law 


<r 2<»V'6| h 2> >0 = <^| l, 2> 

evaluated at r 2 . 

The relations analogous to equations (35) to (37) are 


lh_\ = — 

v 2 ~ v w ! 7rP re( r 2)/ 1 I 4Z \ 

1 2 ' p r > 0 77 

v 4p y Trr 2 j 


<'V‘«| h 2> = ^ 


and 


v w ~ v 2 _ ^2 P r0l 

[ r i)i 


v 4pZ \ 

\ 'W 


where V 2 is the flow velocity at and 

, ( 2 - ° l) 1 

?2 = 


(43) 


(44) 


(45) 


(46) 


(47) 


Equations (37) and (46) provide the integration constants for equation (23), the solu- 
tion of which can be written 


v = 





+ r 2 (r 2 



(48) 


A velocity profile thus exists which conforms both to the simple first-order distribution 
function of equation (22) and to the complete set of boundary conditions listed in equa- 
tions (28) to (32) at r = r^ and a similar set at r = r hence, the present description 
is valid throughout the gas. In particular, the assumption is justified that microscopic 
and macroscopic relations derived in the midst of the gas apply also in the immediate 
proximity of the cylinder surfaces. 

The remainder of this report is concerned with the computation of the torque exper- 
ienced by the inner cylinder. 
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TORQUE ON THE INNER CYLINDER 


As is evident from equation (1), the element P r g of the pressure tensor is the 
essential parameter for the computation of the torque exerted on the inner cylinder by 
the shearing action of the rotating gas. The results of the preceding section are sufficient 
to determine this parameter through linear terms in the flow velocity and cubic terms in 
the mean free path; in particular, equations (20), (37), (46), and (48) combine to yield 


P r0 " 


8p^rjr 2 


2r l r 2 


( r 2 2 - 'l 2 ) : 


Wre (*i) I 1 + if^_) + VA. (-2) f 1 - if Q - 4P ' riVw 


7 TV 


(49) 


Accordingly, p r g( r i[) and p r 0( r 2) satisf y the simultaneous expressions 


*1*1 t 

\ + «* \ 

2 ?o r 2 / 

p ( r \ | ^ 6 ( 

U 2 \ 

\ p / r» \ = 

8pfr 2 v w 

( r 2 2 ’ r l 2 ) 


v ^i r i/ 

l) 2 r 2 1 

r 2 r l 

\ "*2*2/ 

' r0 ^ 2 ) 7TV j 

( r 2 2 - ^ 2 ) 


(50) 


and 


2 ?l r l i 

1 + ^ ^ 

| p re( r l) + 

— 

1 + - 

2 ^2 r l 2 

1 

_ «* \ 

p /r \ = 

8pfri 2 v w 


r 2 2 - ' 

V ^i r i) 

x *r 

! 

[ r 2 2 - r l 2 ) 

l r 2 

V 7T *2 T 2/ 

re( r 2) 

v ’ 7TV | 

f r 2 _ r 2 N i 
^ r 2 r l ) 

l r 2 


(51) 


Finally, the combination of equations (1) and (17) with the angular frequency 


CjO = 


v w 


(52) 


and the solution of equations (50) and (51) for gives 

L = 47n ^ r l 2r 2 2 L _ 2 (^ l r 2 3 + Vl 3 ) 


- r. 


l r 2 ( 


r„ 2 - r 


2 ) 


4 ( g l r 2 3 * ?2 r l 3 ) 


r 2 r 2 
r l r 2 


fc 2 - n 2 )‘ 


> . - r 2 2 )(r/- ri 2 )‘ 

77 ( ? l r 2 3 + ? 2 r l 3 ) 


(53) 


15 



3 

The justification for retaining l terms follows from the first equality in equation (49) 
and the observation that and Vg are correct through since the leading contri- 
bution to P r g appearing in equations (37) and (46) is proportional to l. 

Besides correcting in a general manner the existing expansion (see ref. 1, p. 298) 


L = 


Airrjcor^r^ 




2 ( g l r 2 3 ^ 2 r l 3 ) 
r l r 2( r 2 2 - r l 2 ) 


1 


9 2 

4 7777 cor ^rg 



2 (b r 2 3 + ?2 r l 3 ) + 4 (? l r 2 3 ^2 r l 3 ) 2 

r l r 2( r 2 2 - r X 2 ) rj2r 2 2 (r 2 2 - rj2) 2 


(54) 


- r ]^, equation (53) is applicable for 
even greater rarefication if the cylinder separation d = rg - r^ is small compared 
with r^. In particular, the neglect of quadratic and higher powers of d/r^ in equa- 
tion (53) yields 


through cubic terms in the Knudsen number l j i 


277 Tj cor 


?l + ?2 



1 + 



(55) 


which corresponds to the drag force per unit area 


D = 



? 1 + ^2 


-1 


(56) 


in plane Couette flow (ref. 3) as r^ approaches infinity. 

Since the approach of r^ to infinity, or more precisely the neglect of r^ com- 
pared with unity, makes equation (22) an exact solution of the velocity-linearized kinetic 
equation, the following reduction of equation (55) is limited only by negligible values of 
Ajr^ and small values of (i /d) : 


L = 


2 77?/ cor ^ 
d 



h 



(57) 


No specific restriction is placed on l/d by itself; hence, for sufficiently small Mach 
numbers of the outer cylinder, equation (57) is valid for Knudsen numbers extending well 
into the transition flow regime. Actually, the limitation on d/r ^ need not be quite as 
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strict as stated because the contribution to L from the term in equation (22) which does 
not survive the approach of r^ to infinity also does not survive the elimination of quad- 
ratic and higher powers of d/r^ in equation (53). Thus equation (55) is no more 
restricted than equation (57) as far as the Knudsen number is concerned and permits 
larger values of d/r^. 

If and are eac h equated to £, equation (57) becomes identical with the 
very early (1913) expression deduced by Timiriazeff (ref. 1, p. 298). The preceding 
discussion explains why this formula was found to hold closely even down to pressures 
where it was expected to fail on the basis of less rigorous kinetic theory. Validity cri- 
teria available at that time were based on the Knudsen number alone instead of the quan- 
tity (Z/d)^v w/ /v^ suggested by equations (20) and (22). 

Another consequence of equation (53) is that low-pressure reductions in the torque 
are enhanced by the same geometric effect discussed in the paragraph following equa- 
tions (37) and (38). For example, 


9 2 

Afnrjwr^ 



8 V 




(58) 


in the absence of the conventional contribution to the velocity slip, that is, with 
and 0^2 equated to 2 or ^ and Cg equated to zero. The physical explanation of this 
effect, which is not predicted by the Grad 13-moment approximation because it derives 
from seldom-used elements of the complete heat-flux tensor, is again not apparent. 

Nevertheless, the geometric effect may be very important at low pressures and for 
the intercylinder separation a sizable fraction of r^. The significance is illustrated by 
the following sample calculations for rg = 2r ^ and the slip distances and ^ 
equal to l: 


L = 


1 67777 cor.. 



(59) 


from equation (54) and 


L = 


1677 ?7 cor H 



(60) 


from equation (53). A decrease of 35 percent occurs in the last term in parentheses. 
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CONCLUDING REMARKS 


The velocity-linearized Boltzmann equation is solved in the present research for 
the problem of viscous gases flowing between two relatively rotating coaxial cylinders. 
Boundary conditions corresponding to the conservation of particles, momentum, and 
energy are used in place of the obsolete Maxwellian hypothesis that the reflected particles 
can be divided into diffuse and specular groups. Unlike the results of Wang Chang and 
Uhlenbeck, who employed the Maxwellian hypothesis and thus restricted the accommodation 
coefficients to a single value for all moments, no nonanalytic (in the mean free path) con- 
tributions to the velocity distribution function are found. Hence, the use of the mean free 
path as a perturbation expansion parameter is valid, at least for the terms retained in 
this report which is the complete specification through linear terms in the Mach number 
and quadratic terms in the Knudsen number. 

Calculated results include the generalization to cubic terms in the Knudsen number 
of the Navier-Stokes expression for the torque exerted on the inner cylinder by the 
shearing action of the flowing gas. Grad's 13-moment velocity distribution function is 
shown to be insufficient for this purpose; in particular, seldom-used elements of the com- 
plete heat-flux tensor are required unless the separation between the cylinders is small 
compared with their radii. If this last condition is satisfied and the relative Mach number 
of the rotating cylinders is sufficiently small, a major restriction on the mean free path 
is removed and the present formulas are valid for Knudsen numbers extending well into 
the transition flow regime. 

Besides defining and comparing the limitations on various expressions for the 
torque, the present research introduces contributions to the velocity slip and the torque 
which, except for a shear factor in the slip, act independently of the gas-surface accom- 
modation coefficients for parallel momentum absorption. Although the physical interpre- 
tation of these contributions is not clear, their effects can be very substantial at low pres- 
sures and for intercylinder separations which are sizable fractions of the radii. They 
are not predicted at all by the 13-moment distribution function and are predicted incor- 
rectly by previous applications of the Burnett equations. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 9, 1970. 
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